环境质量分析的数理统计方法

我们掌握的环境数据都是从环境总体抽样出来的样本，必须有一个科学方法去处理样本才能获得对环境总体的正确认识，正确指导环境决策，才能从样本中获取尽可能多的信息。     

Spatial and space-time scan statistics to detect low rate clusters of sex ratio

Whether general environmental exposures to endocrine disrupting chemicals (including pesticides and dioxin) might induce decreased sex ratios (male/female ratio at birth) is discussed. To address this issue, the authors looked for a space-time clustering test which could detect local areas of significantly low risk, assuming a Bernoulli distribution. As a matter of fact, if the endocrine disruptor hypothesis holds true, and if the sex ratio is a sentinel health event indicative of new reproductive hazards ascribed to environmental factors, then in a given region, either a cluster of low male/female ratio among newborn babies would be expected in the vicinity of polluting municipal solid waste incinerators (MSWIs) (supporting the dioxin hypothesis), or local clusters would be expected in some rural areas where large amounts of pesticides are sprayed. Among cluster detection tests, the spatial scan statistic has been widely used in various applications to scan for areas with high rates, and rarely (if ever) with low rates. Therefore, the goal of this paper was to check the properties of the scan statistics under a given scenario (Bernoulli distribution, search for clusters with low rates) and to assess its added value in addressing the sex ratio issue. This study took place in the Franche-Comté region (France), mainly rural, comprising three main MSWIs, among which only one had high dioxin emissions level in the past. The study population consisted of 192,490 boys and 182,588 girls born during the 1975–1999 period. On the whole, the authors conclude that: (i) spatial and space-time scan statistics provide attractive features to address the sex ratio issue; (ii) sex ratio is not markedly affected across space and does not provide a reliable screening measure for detecting reproductive hazards ascribed to environmental factors.     

Modelling the distribution of plant species using the autologistic regression model

For modeling the distribution of plant species in terms of climate covariates, we consider an autologistic regression model for spatial binary data on a regularly spaced lattice. This model belongs to the class of autologistic models introduced by Besag (1974). Three estimation methods, the coding method, maximum pseudolikelihood method and Markov chain Monte Carlo method are studied and comparedvia simulation and real data examples. As examples, we use the proposed methodology to model the distributions of two plant species in the state of Florida.     

Development of an empirical equation for the transverse dispersion coefficient in natural streams

In this study, a new empirical equation for the transverse dispersion coefficient has been developed based on the hydraulic and geometric parameters in natural streams using a regression technique. First, a total of 32 data sets in 16 streams were collected. Among those sets, 16 sets were used for deriving the new equation, and the other 16 sets were used for verifying the equation. Then, through dimensional analysis, it was found that the normalized transverse dispersion coefficient is associated with several parameters such as sinuosity, aspect ratio, and a friction term. The robust least square method was applied to estimate regression coefficients. The newly proposed equation was proven to be superior in explaining the dispersion characteristics of natural streams more precisely compared to the existing equations.     

Efficient statistical mapping of avian count data

We develop a spatial modeling framework for count data that is efficient to implement in high-dimensional prediction problems. We consider spectral parameterizations for the spatially varying mean of a Poisson model. The spectral parameterization of the spatial process is very computationally efficient, enabling effective estimation and prediction in large problems using Markov chain Monte Carlo techniques. We apply this model to creating avian relative abundance maps from North American Breeding Bird Survey (BBS) data. Variation in the ability of observers to count birds is modeled as spatially independent noise, resulting in over-dispersion relative to the Poisson assumption. This approach represents an improvement over existing approaches used for spatial modeling of BBS data which are either inefficient for continental scale modeling and prediction or fail to accommodate important distributional features of count data thus leading to inaccurate accounting of prediction uncertainty.     

MULTIVARIATE TREND TESTING OF LAKE WATER QUALITY1

ABSTRACT: Multivariate methods of trend analysis offer the potential for higher power in detecting gradual water quality changes as compared to multiple applications of univariate tests. Simulation experiments were used to investigate the power advantages of multivariate methods for both linear model and Mann-Kendall based approaches. The experiments focused on quarterly observations of three water quality variables with no serial correlation and with several different intervariable correlation structures. The multivariate methods were generally more powerful than the univariate methods, offering the greatest advantage in situations where water quality variables were positively correlated with trends in opposing directions. For illustration, both the univariate and multivariate versions of the Mann-Kendall based tests were applied to case study data from several lakes in Maine and New York which have been sampled as part of EPA's long term monitoring study of acid precipitation effects.     

ANTIDEGRADATION WATER QUALITY CRITERIA FOR THE DELAWARE RIVER: A DISTRIBUTION-FREE STATISTICAL APPROACH1

ABSTRACT: Existing water quality for the Middle Delaware Scenic and Recreational River is significantly better than is required by current standards, leaving a potential for degradation. A method is presented for deriving candidate antidegradation water quality criteria for this segment of the Delaware River using statistical analysis of historic (ambient) water quality data. Data for 34 water quality parameters are first evaluated for data density, serial correlation, trend, seasonality, and other factors. These preliminary analyses are based on observation of data plots and application of distribution-free statistical techniques that are insensitive to outliers and are robust to relatively mild violations of basic assumptions. Data for 12 of the parameters have sufficient density for further analysis and can reasonably be modeled as independent and identically distributed over time (either seasonally or for the entire data sets). For these 12 parameters, distribution-free statistical methods are developed and used to derive intervals within which there is high confidence (usually greater than 95 percent) that the quantiles with potential use as anti-degradation criteria (the 0.85^th, 0.90^th, and 0.95^th quantiles in this study) for a particular parameter lie.     

REGIONALIZATION OF LOW-FLOW FREQUENCY ESTIMATES: AN ALABAMA CASE STUI)Y1

ABSTRACT: Low-flow estimates, as determined by probabilistic modeling of observed data sequences, are commonly used to describe certain streamflow characteristics. Unfortunately, however, reliable low-flow estimates can be difficult to come by, particularly for gaging sites with short record lengths. The shortness of records leads to uncertainties not only in the selection of a distribution for modeling purposes but also in the estimates of the parameters of a chosen model. In flood frequency analysis, the common approach to mitigation of some of these problems is through the regionalization of frequency behavior. The same general approach is applied here to the case of low-flow estimation, with the general intent of not only improving low-flow estimates but also illustrating the gains that might be attained in so doing. Data used for this study is that which has been systematically observed at 128 streamflow gaging sites across the State of Alabama. Our conclusions are that the log Pearson Type 3 distribution is a suitable candidate for modeling of Alabama low-flows, and that the shape parameter of that distribution can be estimated on a regional basis. Low-flow estimates based on the regional estimator are compared with estimates based on the use of only at-site estimation techniques.     

IMPROVED ESTIMATION OF ECOLOGICAL EFFECTS USING AN EMPIRICAL BAYES METHOD1

ABSTRACT: The US Environmental Protection Agency is currently developing guidance for assessing environmental impact using biocriteria within a regulatory framework. Of particular interest are statistical methods of design and analysis to test for impairment of biological assemblages in stream ecosystems associated with water pollution. Current recommendations emphasize frequentist statistical methods with the problem expressed as one of classical hypothesis testing. An empirical Bayes approach is proposed here as an alternative and applied for multi-site inference. The advantages of an empirical Bayes approach, in particular the substantive contribution of collateral information from nearby sites, are discussed in contrast to traditional methods that employ site-specific information only. The approach is illustrated in an application concerning trends in the Index of Biotic Integrity (IHI) for the Scioto River in Ohio.     

A WATERSHED-LEVEL ECOLOGICAL RISK ASSESSMENT METHODOLOGY1

ABSTRACT: We present an ecological risk assessment methodology at the watershed level for freshwater ecosystems. The major component is a pollutant transport and fate model (a modified EUTROMOD) with an integrated uncertainty analysis utilizing a two-phase Monte Carlo procedure. The uncertainty analysis methodology distinguishes between knowledge uncertainty and stochastic variability. The model assesses the ecological risk of lentic (lake) ecosystems in response to the stress of excess phosphorus resulting in eutrophication. The methodology and model were tested on the Wister Lake watershed in Oklahoma with the lake and its trophic state as the endpoint for ecological risk assessment. A geographic information system was used to store, manage, and manipulate spatially referenced data for model input.